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基于 Edgeworth 级 数 和 多 重 线 性 化 的 
概率 潮流 算法 
( 广西 电网 有 限 责任 公司 电力 科学 研究 院 ” ”南宁 530023 ) 


摘要 ， 对 含 随机 变量 的 大 规模 概率 潮流 问题 ， 若 使 用 Gram-Charlier 级 数 展开 式 逼 
近 随 机 变量 的 概率 分 布 ， 计 算 量 较 大 。 建 立 考虑 风电 出 力 随机 性 和 负荷 波动 性 的 概率 
潮流 模型 ， 并 采用 Edgeworth 通过 的 Hermite 正 交 和 矩阵 的 递 推 性 降低 级 数 高 阶 展开 的 
计算 复杂 度 ， 针 对 级 数 变换 法 高 次 项 被 舍 去 影响 有 逼近 精度 的 问题 ， 采 用 多 重 线性 化 的 
方法 通过 对 系统 总 有 功 功率 的 均匀 划分 ， 减 小 因 谢 流 方程 线性 化 时 输入 的 随机 变量 的 
变化 范围 较 大 而 引起 的 截断 误差 。IEEE 39 布点 系统 的 算 例 仿真 结果 验证 了 本 方法 的 
有 效 性 。 

关键 词 ， 概率 滑 流 ”多 重 线性 化 ” 递 推 性 ”Edgeworth 级 数 
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A Probabilistic Power Flow Algorithm Based on Edgeworth 
Series and Multiple Linearization 


Lianeg .ie 
( Electric Power Research Institute of Guangxi Power Grid Co.,Ltd. 
Nanning 530023 China ) 


Abstract: On the problem of large-scale probabilistic power flow problems with 
random variables, if the Gram-Charlier series expansion was used to approximate the 
cumulative distribution of random variables, the amount of calculation was large. In order 
to reduce the computational complexity of higher order expansion by using the recurrence 
of Hermite orthogonal matrix of edgeworth, a probabilistic power flow model considering 
the randomness of wind power output and load fluctuation is established. According to the 
method of high-order series transform is eliminated the influence of calculation accuracy, 
the method of multiple linear systems by uniform division of the total active power is 
used to reduce the truncation error caused by the large variation range of the input random 
variables due to linearization of the power flow equation. Simulation results of IEEE 39 
system verify the effectiveness of the proposed method. 
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概率 潮流 分 析 可 以 帮助 调度 运 维 人 员 评估 系 统 
随机 因素 变化 对 潮流 安全 稳定 运行 的 影响 ， 以 确定 
安全 调度 策略 。 在 概率 潮流 分 析 中 ， 除 了 传统 火电 
机 组 的 出 力 和 网 络 阻抗 等 可 预测 因素 外 ， 还 需 考 虑 
风电 等 新 能 源 的 不 稳定 出 力 和 负荷 波动 等 随机 扰动 
因素 对 系统 安全 性 的 影响 中。 

文献 [2] 根据 节点 电压 和 支 路 潮流 的 期 望 值 及 
灵敏 度 矩 了 泗 ， 计 算 了 负荷 、 风 电机 组 出 力 等 的 半 
不 变量 ， 由 Gram-Charlier 级 数 展开 求 得 概率 分 布 
函数 ， 但 Gram-Charlier 展开 需要 通过 各 模拟 抽样 
序列 计算 展开 式 的 各 阶 距 ， 在 高 阶 展开 时 计算 量 
较 大 。 文 献 [3] 采 用 蒙特 卡 洛 仿真 (Monte Carlo 
Simulation，MCS) 通过 简单 随机 采样 得 到 输入 变 
量 的 样本 ， 再 对 每 个 采样 点 进行 确定 性 潮流 计算 ， 
最 后 统计 获得 各 状态 量 的 分 布 情况 ， 该 方法 原理 简 
单 ， 在 保证 样本 规模 足够 大 的 情况 下 能 够 获得 很 高 
的 精度 ， 缺 点 是 仿真 次 数 多 ， 在 用 于 求解 大 规模 系 
统 的 概率 潮流 时 耗 时 较 长 。 

本 文 首先 建立 考虑 风电 出 力 不 稳 定性 和 负荷 波 
动 的 交流 概率 潮流 模型 ， 并 使 用 多 重 线性 化 的 方法 
通过 对 系统 总 有 功 功 率 的 均匀 划分 ， 将 潮流 方程 线 
性 化 ， 然 后 通过 Edreworth 级 数 展 开 计 算出 潮流 方 
程 中 随机 状态 变量 的 概率 分 布 。 最 后 通过 IEEE 39 
节点 系统 的 算 例 仿 真 分 析 了 风电 出 力 的 不 确定 性 和 
负荷 波动 对 基于 级 数 展开 的 概率 潮流 逼近 精度 的 影 
响 ， 并 比较 了 不 同 级 数 展开 的 逼近 精度 和 计算 速度 。 


2 基于 Edgeworth 变换 和 多 重 线性 化 的 
概率 潮流 


2.1 交流 潮流 模型 及 其 半 不 变量 


VOYV,cosd,+P,-P,-P,,=0 iESwNSo, 
jEi 


沁 
1 (1) 
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式 中 ，Spv 和 Spo 分 别 为 系统 PV 和 PQ 节点 集合 ，; 
Po 和 0 为 节点 i 火电 机 组 发 出 的 有 功 、 无 功 功 
率 ; Pw 和 Ow; 为 节点 i 处 风电 机 组 输出 的 有 功 、 无 
功 功 率 ，Pi; 和 01; 为 节点 i 负荷 的 有 功 、 无 功 功率 ; 
太 和 656; 为 节点 i 电压 幅 值 和 相 角 ;为 节点 导 纳 条 
阵 元 素 ; aj 为 相应 的 相 角 ，6, =6;-6j-ay。 本 文 
概率 潮流 萎 虑 负荷 的 波动 性 与 风电 机 组 出 力 的 不 确 


ChinaXiv 合 作 期 刊 


定性 ， 则 Pw;、Qw、Pi、Q4 为 随机 变量 。 

将 式 (1) 的 放流 方程 线性 化 ， 假 设 各 节点 注入 
功率 的 随机 变化 是 独立 的 ， 则 系统 状态 变量 实际 上 
是 各 注入 功率 随机 变量 的 线性 和 ， 且 权重 系数 为 灵 
敏 度 系数 "。 基 于 此 ， 将 式 (1) 在 基准 运行 点 处 采 
用 泰勒 展开 并 忽略 二 阶 以 上 的 高 次 项 ， 可 得 其 单 点 
展开 模型 为 


J 02) 


Z=Z,+T,W-W) 


式 中 ， 各、Zo、 克 分 别 为 基准 运行 点 处 的 市 点 电压 、 
支 路 功率 和 节点 注入 功率 ， 到 为 市 点 注入 功率 的 随 
机 变量 ，S6 与 Te 分 别 为 节点 电压 和 支 路 功率 对 注入 
功率 变化 的 灵敏 度 和 矩阵 ， 且 So= ji ，7m= Go ， 
其 中 万 为 雅 可 比 矩 阵 ，GCo=(9ZwWaXo)lx-x。 
运用 半 不 变量 的 特性 可 以 得 到 状态 变量 AX 和 
AZ 支 路 潮流 的 各 阶 半 不 变量 ， 即 ” 
J 


AZ® 全 T® x AW® (3) 
0 


式 中 ，So® 和 To 分 别 为 矩阵 5S。 和 7 中 元 素 的 次 
曙 构 成 的 矩阵 ，AW 为 节点 注入 功率 变化 量 的 次 
需 ， 从 而 可 以 根据 各 节点 发 电机 出 力 、 风 机 出 力 和 
负荷 分 布 的 情况 下 求 出 发 电机 与 负荷 的 半 不 变量 。 
2.2 概率 潮流 方程 的 Edgeworth 级 数 展开 

传统 的 Gram-Charlier 级 数 展 开 用 于 逼近 待 求 
状态 变量 和 支 路 潮流 的 累积 分 布 函 数 (Cumulative 
Distribution Function，CDF) 时 ， 在 每 次 迭代 时 需 
要 逐 项 计算 各 阶 Hermite 正 交 矩阵， 计算 量 较 大 。 
对 此 ， 采 用 Edgeworth 展开 法 ， 可 以 通过 标准 正 态 
分 布 国 数 逼 近 服 从 任意 分 布 的 随机 变量 的 概率 分 布 
函数 ， 且 渐进 展开 的 方式 计算 较 简 便 包 。 首 先 定义 
随机 变量 x 的 标准 形式 为 ?， 则 

ye (4) 

式 中 , 4 为 x 的 均值 ,，o 为 标准 差 ， 设 F(x) 和 fx) 
分 别 为 的 CDF 和 概率 密度 函数 (Probability 
Density Function, PDF)， 有 是 fx)=(x)， 则 x) 的 
Edgeworth 展开 表示 为 


f(D = p07) OD) HD)+ ee HO)+ 
k kk 
+ “H,(y) + (5) 
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式 中 ,gp0y) 为 标准 正 态 分 布 的 PDF，Hb) 为 i 阶 
Hermite 正 交 矩阵, ;为 y 的 i 阶 半 不 六 量 ， 且 . 


hE -eo0)] (6) 
式 中 ，E(: ) 表示 期 望 值 ，H,Qy) 有 如 下 递 推 关系 


0 而 
H(y)=1,H(y)=» 

由 于 Edgeworth 的 Hermite 矩阵 在 Hilber 空间 
上 具有 正 交 性 ， 由 式 (7), 的 后 项 可 由 前 项 的 递 
推 关 系 确定 下 ， 求 出 前 项 后 ， 后 项 可 根据 前 项 的 多 
项 式 技术 快速 获得 ， 从 而 减少 计算 量 。 
2.3 潮流 方程 的 多 重 线性 化 

但 当 输 入 的 随机 变量 的 变化 范围 较 大 时 ， 由 于 
级 数 展开 后 的 高 次 项 被 舍 去 ， 采用 式 (2) 的 单 点 线 
性 化 模型 会 引起 较 大 的 截断 误差 ， 为 了 在 提高 级 数 
变换 法 计算 速度 的 同时 不 至 于 损失 太 多 精度 ， 下 面 
引入 一 种 多 重 线性 化 的 方法 。 设 Pi 为 系统 总 的 有 
功 功率 ， 有 


P=>P->P (8) 
之 


式 中 ，NN 为 PQ 节点 个 数 ， 己 为 节点 站 负荷 有 功 功 
率 ，Na 为 PV 节点 个 数 ，P 为 发 电机 节点 /的 有 功 
输出 。 

由 于 负荷 功率 和 新 能 源 出 力 均 为 随机 变量 ， 可 
知 Po 也 是 随机 变量 。 多 重 线性 化 的 方法 首先 将 Ps 
等 间距 地 划分 为 Ri 个 区 域 。 在 区 域 R ~ Ri 内 选取 
相应 的 基准 运行 点 ， 并 将 潮流 方程 在 各 点 处 分 别 线 
性 化 ， 即 


Xi = + Sr (Wi -Wi,) 
Xi Xo0 + (Wi —Wi,) 


Xr = Xro + Sr (Wher 9) 


RT 


式 中 ，Si; 为 第 Ri 个 区 域 对 应 的 灵敏 度 矩 阵 ，Xi 为 
第 Ri 个 区 域 的 状态 向 量 ， 克 ;为 第 Ri 个 区 域 的 基准 
功率 向 量 。 

综 上 ， 本 文 基于 多 重 线性 化 和 Edgeworth 级 数 
展开 的 概率 潮流 算法 流程 如 图 1 所 示 。 


3 ” 算 例 分 析 
为 验证 本 文 所 提 方 法 的 可 行 性 ， 以 IEEE 39 节 
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输入 基础 数据 


本 
计算 出 节点 电压 和 支 路 潮流 的 期 望 值 
加 ，Zo 及 灵敏 度 和 矩阵 $。，T 


将 系统 总 有 功 功率 等 间距 地 划分 ， 在 各 
分 区 的 基准 运行 点 将 户 流 方程 线性 化 


L 1 
根据 分 布 率 求 负荷 的 各 阶 | | 根据 分 布 率 求 风电 机 组 的 
半 不 变量 各 阶 半 不 变量 


介 半 不 变量 


了 
求 各 节点 注入 功率 变化 量 
AW 的 各 阶 半 不 变量 


浊 
分 别 计算 状态 变量 AX 和 支 
路 潮流 AZ 的 各 阶 半 不 变量 


下 
1Edgeworth 级 数 展开 求 得 
AX 和 AZ 的 PDF 和 CDF 


1 算法 流程 
Fig.1 Flow chart ofthe proposed method 


点 系统 为 算 例 ， 该 系统 含 39 个 节点 ，46 条 支 路 ， 假 
设 风机 出 力 服 从 Weibull 分 布 ， 负 荷 服 从 正 态 分 布 ， 
各 随机 变量 之 间 相 互 独立 ， 模 型 见 文献 [10]， 且 以 
原 系 统 负 荷 预 测 数 据 为 分 布 期 望 值 。 

为 验证 本 文 方法 的 准确 性 ,分 别 将 本 文 方法 和 
Gram-Charlier 级 数 概率 潮流 计算 结果 与 基于 MCS 
的 结果 进行 比较 ，MCS 法 抽样 次 数 Ni 取 1 000 次 ， 
采用 方差 和 根 均 值 (Average Root Mean Square， 
ARMS)"" 作为 逼近 精度 评价 指标 ， 即 


N 
(Ses -CMC,) 
ARMS=A (10) 
N 


式 中 ，CCE 和 CMC, 分 别 为 半 不 变量 法 和 MCS 法 
得 到 的 输出 随机 变量 累积 分 布 曲线 上 第 i 个 点 的 值 ， 
NN 为 市 点 数 。 

分 别 对 两 个 场景 进行 分 析 : 场景 1 只 考虑 负 和 蓓 
波动 ， 场 景 2 考虑 负 符 波动 和 风电 机 组 出 力 的 波动 。 

图 2 显示 了 场景 1 的 节点 电压 幅 值 精度 指标 的 
均值 和 最 大 值 随 负 和 荷 波动 标准 差 的 变化 情况 。 可 见 ， 
当 系 统 只 考虑 负荷 波动 时 ， 两 种 级 数 展开 的 半 不 变 
量 法 台 近 精度 均 较 低 ， 这 是 由 于 本 文 假设 负荷 波动 
服从 正 态 分 布 ， 所 述 两 种 级 数 都 是 将 输入 随机 变量 
的 分 布 函数 表示 为 由 正 态 随机 变量 各 阶 导 数组 成 的 
级 数 '"， 故 两 者 逼近 精度 均 较 高 ， 当 负荷 波动 性 加 
大 后 ， 半 不 变量 法 逼近 精度 减 小 。 这 是 由 于 负荷 波 
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一 -Gram-Charlier 级 数 
--e。- Edgeworth 级 数 
最 大 值 


ARMS(%) 
SS 
ULD 


均值 
| 10 13 20 25 30 
负 和 蓓 波动 标准 差 (%) 
图 2 场景 1 节点 电压 幅 值 的 ARMS 均值 和 最 大 值 
Fig.2 Average and maximal ARMS of bus voltage amplitude 


incasel 

动 标 准 差 较 大 时 ， 加 大 了 注入 功率 远离 负荷 功率 期 
望 值 的 概率 ， 使 得 潮流 线性 化 处 理 引 起 的 误差 随 之 
增 大 。 

场景 2 在 场景 1 的 基础 上 ， 将 Vestas V112 一 
2.5MW 型 风电 机 组 接 入 市 点 11， 切 入 、 切 出 和 和 额 
定 风 速 分 别 为 4、25 和 16m/s， 采 用 恒 功 率 因 数控 
制 策略 。 风 速 的 双 参 数 Weibull 分 布 参 数 分 别 取 = 
3.97 和 c=10.7。 

场景 2 两 种 方法 在 不 同 风电 机 组 容量 下 的 电压 
幅 值 逼近 精度 指标 如 图 3 所 示 。 由 图 3 可 见 : @ 当 
风电 机 组 接 入 容量 增加 时 ， 两 种 方法 的 ARMS 指标 
值 均 增 大 ， 即 半 不 变量 法 逼近 精度 呈 降 低 趋势 ，@ 
由 于 风电 机 组 容量 对 应 为 随机 变量 取 值 的 最 大 值 的 
范围 ， 故 ARMS 最 大 值 随 风电 机 组 容量 增 大 的 趋 
势 较 明显 ， 而 均值 为 随机 变量 每 个 采样 点 根据 采样 
次 数 平均 化 的 结果 ， 由 于 采样 次 数 基 数 较 大 ， 故 均 
值 变 化 的 趋势 较 平 缓 ，@ 本 文 基 于 多 重 线性 化 的 
Edgeworth 级 数 展 开 法 的 拟 合 精度 比 Gram-Charlier 
级 数 稍 好 ， 但 区 别 不 大 ， 当 风电 机 组 装机 容量 增 大 
时 ， 区 别 逐 浙 增 大 。 这 是 由 于 本 文采 用 多 重 线 性 化 
方法 ， 将 系统 总 的 有 功 功 率 等 间距 地 划分 为 若干 个 
区 域 ， 在 各 区 域 选取 基准 运行 点 将 潮流 方程 线性 化 ， 
可 减少 采用 单 点 线性 化 模型 使 得 基准 点 和 实际 运行 
点 之 间 偏 差 导 致 的 截断 误差 。 


6 


一 Gram-Charlier 级 数 
--e- Edgeworth 级 数 


最 大 值 


ARMS(%) 
LD 


0 5 10 15 20 25 
风电 机 组 容量 /MW 


图 3 场景 2 节点 电压 幅 值 的 ARMS 值 和 最 大 值 
Fig.3 Average and maximal ARMS of bus voltage amplitude 


in case 2 
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为 比较 本 文 方法 、Gram-Charlier 展开 法 与 
MCS 方法 的 求解 耗 时 ， 在 场景 2 的 风机 接 入 容量 为 
25MW 时 ， 计 算 线 路 13-14 上 的 有 功 潮流 概率 分 布 ， 
结果 见 下 表 ，MCS 法 不 作 自 我 比较 。 可 见 ， 本 文 方 
法 的 逼近 精度 稍 低 于 Gram-Charlier 展开 法 ， 这 是 由 
于 本 文 的 Edgeworth 展开 法 仅 通过 Hermite 矩阵 信 
息 近 似 拟 合 高 阶 项 ， 存 在 拟 合 误差 。 当 采样 次 数 为 
1 000 时 ， 本 方法 计算 速度 比 MCS 法 快 约 22%， 这 
是 由 于 MCS 法 要 生成 相应 分 布 类 型 的 随机 样本 反 
复 进行 模拟 场景 的 确定 性 潮流 计算 ， 故 计算 量 较 大 ， 
本 方法 计算 速度 比 Gram-Charlier 展开 法 稍 快 ， 这 是 
由 于 Gram-Charlier 展开 法 每 阶 展开 式 均 要 重新 计算 
各 阶 距 ， 而 Edgeworth 多 项 式 在 Hilber 空间 上 具有 
正 交 性 ， 它 的 高 阶 Hermite 矩阵 可 由 低 阶 项 的 递 推 
关系 确定 ， 求 出 低 阶 距 后 ， 高 阶 距 可 根据 低 阶 距 信 
息 快 速递 归 获 得 ， 故 计算 速度 较 快 。 

表 不 同方 法 的 求解 耗 时 对 比 


Tab. Comparison of calculation time of different methods 


本 文 方法 ”Gram-Charlier 级 数 MCS 
ARMS(%) 1.34 1.3 / 
计算 耗 时 As 14.04 17.87 64.54 


4 ”结束 语 


本 文 提出 一 种 基于 Edgeworth 级 数 展 开 的 概率 
潮流 算法 ， 利 用 Edgeworth 多 项 式 各 阶 Hermite 算 
阵 的 递 推 性 降低 高 阶 展 开 的 计算 复杂 度 。 此 外 ， 引 
入 一 种 多 重 线性 化 的 方法 将 潮流 方程 线性 化 ， 有 利 
于 减少 级 数 展开 法 求解 概率 潮流 模型 时 舍 去 高 阶 
项 造成 的 舍 入 误差 。IEEE 39 节点 的 数值 仿真 中 以 
MCS 法 计算 结果 为 参考 ， 以 ARMS 为 评价 指标 ， 
对 比分 析 了 负荷 波动 和 风机 输出 这 两 种 输入 随机 变 
量 对 Edgeworth 以 及 Gram-Charlier 级 数 展 开 下 半 不 
变量 法 概率 潮流 计算 结果 的 准确 性 的 影响 。 结 果 表 
明 ， 风 电 接 入 量 增加 后 半 不 变量 法 的 拟 合 准确 度 均 
减少 ， 但 本 文 方法 的 精度 受 非 正 态 分 布 随机 变量 规 
模 的 影响 相 比 Gram-Charlier 级 数 展开 法 较 小 。 比 较 
了 本 文 方法 与 Cornish-Fisher 展开 和 MCS 法 计算 速 
度 上 的 差别 ， 结 果 表 明 在 ARMS 指标 相差 不 大 时 ， 
本 文 方法 比 Gram-Charlier 展开 法 稍 快 ， 比 MCS 法 
快 约 22%， 体 现 了 本 文 方法 在 计算 速度 上 的 优势 。 
如 何在 保证 概率 潮流 分 布 的 将 近 精 度 的 同时 提高 计 
算 速 度 仍 值得 进一步 研究 。 
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